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We study recent Deeply Virtual Compton scattering (DVCS) data within a dual parameterization 
of the Generalized Parton Distributions (GPDs). This parameterization allows to quantify the 
maximum amount of information, that can be extracted from DVCS data, in a "quintessence" 
function. We present a "zero step" model for the latter solely based on the forward quark density, 
providing a parameter free prediction for the imaginary part of the DVCS amplitude. It is shown 
that the bulk effect of the ep — > ep7 beam helicity cross section difference can be understood within 
such a model. 
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The generalized parton distributions (GPDs) [l[ (see 
Refs. [10, 110] for recent reviews) describe the response 
of the target hadron to the well-defined QCD quark-gluon 
operator of the type : V> Q (0) Pe lg fo dx » Af * ipp(z), which 
is defined on the light-cone, i.e. z 1 = 0. The hard exclu- 
sive processes provide us with a set of new fundamental 
probes of the hadronic structure. The hadron image seen 
by such non-local probes is encoded in the dependence 
of GPDs on its variables. The central quantities that 
contain invaluable information on the nucleon structure 
are momentum transfer (t) dependent quark densities. 
These densities can be obtained as the £ — ► limit of the 
nucleon GPDs H(x, £, t) and E(x, £, t) : 



q(x, t) = lim H(x, £, t), e(x,t) 



KmE(x,Z,t), (1) 



where x ± £ correspond to the quark longitudinal mo- 
mentum fractions. The first density q(x,t) (at t = Q it 
reduces to the usual quark distribution measured in DIS) 
is related to the distribution of quarks and anti-quarks 
in the longitudinal momentum and transverse plane of 
the nucleon, thus providing the 3D image of quarks and 
anti-quarks in the nucleon [6|. The new function e(x,t) is 
crucial for extraction of the angular momentum carried 
by quarks in the nucleon 0]. 

Although the quark densities q(x, t) and e(x, t) are de- 
fined as the simple limit of GPDs one can not perform 
this limit measuring observables for hard exclusive pro- 
cesses. The reason is that the leading order amplitude of 
hard exclusive reactions (we restrict ourselves to DVCS 
and discuss only the GPD H(x,£,t), the discussion for 
E(x, £, t) is analogous) is expressed as : 



J dxH(x,£,t) 
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which the dependence of the GPDs on the variable x 
is "integrated out" . One cannot completely restore the 
GPD H(x,£,t) from ([2]), hence one is not able to per- 
form a "complete imaging" of the target hadron from 
the knowledge of the amplitude and cannot perform the 
limit ([l} to obtain the key quark densities q(x, t) and 
e(x, t). 

The aim of this Letter is to estimate the contribution 
of the quark densities q(x,t) to DVCS observables. This 
is to be considered as a "zero step" to extract the GPD 
information from the data, as it allows us to quantify the 
deviations in terms of genuine non-forward parts of the 
GPDs. For our analysis we employ the dual parameteri- 
zation of GPDs suggested in Ref. Q. 

The dual parameterization is based on a representation 
of parton distributions as an infinite series of t-channel 
exchanges [ll|, allowing to express the GPD H in terms 
of a set of functions Q2u{x,t) (y = 0,1,2,...), for de- 
tails see [8[- We call the functions Q2v{x,t) forward-like 
because at leading order (LO), the scale dependence of 
the functions Q2u{x,t) is given by the standard DGLAP 
evolution equation, so that these functions behave as 
usual parton distributions under QCD evolution. Fur- 
thermore, the function Qo(x,t) is related to the for- 
ward i-dependent quark densities q(x, t) (which reduces 
at t = to the parton densities q(x) measured in DIS) : 



)o{x,t) = [q + q] (x,t) 



dz 

— [q + q](z,t). (3) 



The amplitude is given by the convolution integral in 



The expansion of the GPD H(x,(,,t) around £ = with 
x fixed to the order ^ 2v involves only a finite number of 
functions Q2^{x, t) with /i < v. One can then express the 
amplitude ^ in terms of forward-like functions as [1, 0, 
E3: 



2 



Im A(£,i) 
Re A(£, t) 
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= 2 



iV(a;,t) 



1 _ 2p + x 2 



2£ i T 2 VI +CC 2 



^ / „— JV(s,t) 



1 + 2p + x 2 VI + a; 2 



•4£>(t). 



(5) 



Here £)(t) is the D-form factor (FF) 



D(t) 



n=l 



d n (t) 



dz 



D(z,t) 



(6) 



where D(z,t) is the D-term [12J. One can check [9j that 
the amplitude given by Eqs. (|4I5[) automatically satisfies 
a dispersion relation with the subtraction consta nt g iven 
by the £>-FF, as it should be on general grounds [13l.ll4|. 
The function N(x,t) in Eqs. (|4I5[) is defined as : 



N(x,t) 



(7) 



The information contained in the LO amplitude is in one- 
to-one correspondence with the knowledge of the func- 
tion N(x,t) and the D-FF D(t), because Eq. (gj) can be 
inverted [9(, i.e. the function N can be expressed un- 
ambiguously in terms of the amplitude. This inversion 
corresponds to the Abel transform tomography [l5|. and 
the corresponding inversion equation has the form [9( : 



N(x,t) 



1 x(l-x 2 ) 

~ (J + 2-2)3/2 



d^_ 

£3/2 



2x 
1+x 2 



<ilmA(^)-^IrnA(£,i)}.(8) 

This equation implies that the function N(x,t) con- 
tains the maximal information about GPDs that is pos- 
sible to obtain from the observables. Therefore this func- 
tion can be called the GPD-quintessence function . An- 
other important feature of the expressions (|4I5|) for the 
amplitude is that one can easily single out the contribu- 
tions to the amplitude coming from the forward parton 
densities. Indeed, the first term in the sum JJJ) is given 
by the function Qq which is related to the (^-dependent) 
parton densities by Eq. ([3]). The big advantage of the 



1 See Refs. 0, [Toj | for the detailed discussion of properties and 
physics interpretation of the GPD-quintessence function. 



dual parameterization is that one can clearly separate 
the contribution of the (t-dependent) parton densities 
from genuine non-forward effects encoded in the func- 
tions Q^iQii ■ ■ ■ ■ The authors of Ref. [l6| developed an 
approach that also allows to separate the contribution 
of forward quark densities to observables. Calculations 
of DVCS observables in the dual parameterization were 
presented in however in this paper a "predefined" 
model for forward-like functions Q2 V has been used and 
the "zero step" separation was not discussed. 

In the following we study the separation between 
effects of forward quark densities and genuine non- 
forward contributions. In order to make this separa- 
tion ("zero step") as clean as possible we choose to an- 
alyze recent JLab/Hall A data [18| on the beam helicity 
(in) dependent cross sections in the e+p — > e' +p + j pro- 
cess, as well as recent beam spin asymmetry data mea- 
sured by the CLAS collaboration [20(. We make such 
choice because the beam helicity dependent cross sec- 
tions and beam asymmetries are directly proportional to 
the imaginary part of the DVCS amplitude. This, in 
principle, gives the possibility to apply directly the Abel 
tomography formula ([H]), and requires a measurement at 
several values of xb- Furthermore, we choose the data 
with Q 2 > 2 GeV 2 and -t < 0.3 GeV 2 . In this kinc- 
matical region one can safely neglect contributions of the 
nucleon GPDs E, H and E relative to the contribution of 
H, also one keeps the contribution of higher twists neg- 
ligible. This allows us to make direct conclusions about 
contribution of the quark densities q(x,t). The region of 
small t makes our results insensitive to the modelling of 
the t-dependence of the quark densities. 

In order to make this "zero step" comparison of the 
data with the contribution of forward quark densities we 
have to make assumptions about the i-dependence of the 
quark densities q(x,t). As the model for this dependence 
we use the Regge motivated Ansatz : q(x, t) — q{x) ■ 
x~ a , where q(x) the forward quark distribution, and 
the slope parameter a' = 1.105 GeV -2 was fixed from 
the form factor sum rule in Ref. [l|| . 

In Fig. [TJ we compare our predictions for the polar- 
ized cross section difference for different beam helici- 
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FIG. 1: Azimuthal angular dependence of the e~p — > e~p7 
polarized cross section difference (d 4 a+ — d 4 cr_)/2 for different 
beam helicities. The black bands are JLab/Hall A data 
Dashed (blue) curves : double distribution parameterization 
with profile parameter b — 1 ; solid (red) curves : dual param- 
eterization including only the forward function Qq. 



ties for which the first data have been published by the 
JLab/Hall A Coll. [13]. This observable is directly pro- 
portional to the imaginary part of the DVCS amplitude. 
The data are compared with predictions computed in 
the dual parameterization and in a double distribution 
model. For both models we use the same Regge-type t- 
dependence for forward ^-dependent quark densities. For 
the double distribution we use the profile function with 
parameter b = 1 as in estimates performed in Refs. [3.l2lj|. 
It is seen that the double distribution model yields po- 
larized cross sections that have a tendency to be some- 
what larger than the data. When using only the forward 
function Qq in the dual parameterization, it is seen from 
Eq. ([3]) that the imaginary part of the DVCS amplitude 
yields a parameter free prediction at t = 2 . It is seen 
that this parameter free prediction yields an amazingly 
consistent description of the polarized cross sections. 

The JLab/Hall A Coll. also published results for un- 
polarized cross sections that are shown in Fig. [2] When 



2 We checked that in the t range shown, the dependence on a' is 
much smaller than the difference between the curves. 



FIG. 2: Azimuthal angular dependence of the e~p — > e~p7 
unpolarized cross section. The black bands are JLab/Hall A 
data The dashed blue curves show the double distri- 

bution result with profile parameter 6=1. The results in 
the dual parameterization including only Qo are shown for 
D(0) = (thick solid red curves), D(0) = -4/3 (dashed- 
dotted curves), and -D(O) = +4/3 (dotted curves). The lower 
thin solid green curves are the Bethe-Heitler cross section. 

comparing the unpolarized cross sections with the model 
independent Bcthc-Hcitler result, one sees that the latter 
dominates the cross section at $ = 0, where it makes up 
for about 85 % of the cross section, at —t — 0.23 GeV 2 . 
However, at $ = 180 deg it is more than a factor 2 below 
the data. Although both double distribution and dual pa- 
rameterization models can explain part of the difference 
with the data, no single model is able to simultaneously 
describe the cross section at $ = and $ = 180 deg, 
in line with the finding of [22J . In particular, within the 
dual parameterization in twist-2 approximation, adjust- 
ing the one subtraction constant D{t) does not allow to 
describe this $ dependence, as is also shown on Fig. [5] 
It was also checked that when adding an estimate for the 
non-forward function Qi no agreement can be found ei- 
ther over the whole $ range. Although the region around 
$ = 180 deg yields zero for the beam helicity asymmetry, 
it is very worthwhile to cross-check this puzzle further. 

In Fig. [31 we show recent exclusive measurements 
of e~p — > e~pj beam helicity asymmetries from 
JLab/CLAS [H and JLab/HallA 0. We note that 
for two middle bins in Fig. [31 where both experiments 
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FIG. 3: Azimuthal angular dependence of the e~p — ■» e~p7 
beam helicity asymmetry (d 4 <r+ — d 4 <T_)/(<i 4 <7+ + d 4 cr-) for 
different kinematics. Black bands in the two middle panels 
are data points from JLab/Hall A [18| • Solid circles are data 
points from JLab/CLAS [201] . Curve conventions as in Fig. [2] 



have overlapping kinematics, the data from both experi- 
ments are consistent with each other. When comparing 
both parameterizations with the data in Fig. [3l we note 
that the double distribution model (dashed curves) yields 
asymmetries that lie systematically above the data. The 
dual parameterization model based on the forward func- 
tion Qq (solid curves) yields a good first ("zero step") 
description of the data, given that no parameter was ad- 
justed here. There is a slight tendency for the asymme- 
tries to be overestimated within the dual parameteriza- 
tion which is merely a reflection of the underestimate of 
the unpolarized cross sections as seen in Fig. [2] 

The success of our "zero step" exercise shows that 
the contribution of the forward quark densities q(x, t) 
at small t to the e~p — > e~p7 polarized cross section dif- 
ference yields the bulk effect. Deviations of this predic- 
tion from the data can be fitted by introducing forward- 
like functions Q2, Qi, ■ ■ ■ which describe the genuine non- 
forward effects in GPDs. The dominance of the forward 
quark densities in the imaginary part of the DVCS am- 
plitude is an important observation. It implies that pre- 
cise measurements of DVCS observables in a wider kine- 
matical region would allow us to extract the £-dependent 
quark densities q(x,t) and e(x,t). The latter can be ac- 
cessed by measurements on the neutron. 

In summary, we studied recently obtained high pre- 
cision DVCS data in the valence region within a dual 
parameterization of the GPDs. This parameterization 



allows to extract a quintessence function which contains 
the maximal amount of information that can be extracted 
from DVCS observables. We established a "zero step" 
model for this function solely based on the forward quark 
density, which provides a parameter free prediction for 
the imaginary part of the DVCS amplitude. It was shown 
that the bulk effect of the ep — > ep-f beam helicity cross 
section difference can be understood within such a model. 
By systematically studying deviations between the data 
and such model will allow to reveal the non-forward ef- 
fects to this quintessence function. 
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